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Abstract. In this paper a simple model for the evolution of the forward 
density of the future value of an asset is proposed. The model allows for 
a straightforward initial calibration to option prices and has dynamics that 
are consistent with empirical findings from option price data. The model is 
constructed with the aim of being both simple and realistic, and avoid the 
need for frequent re-calibration. The model prices of n options and a forward 
contract are expressed as time-varying functions of an (n -|- l)-dimensional 
Brownian motion and it is investigated how the Brownian trajectory can be 
determined from the trajectories of the price processes. An approach based 
on particle filtering is presented for determining the location of the driving 
Brownian motion from option prices observed in discrete time. A simulation 
study and an empirical study of call options on the S&P 500 index illustrates 
that the model provides a good fit to option price data. 



1. Introduction 

Consider a financial market consisting of a collection of European options with 
maturity T > 0, written on the value St of an asset at time T. Suppose that the 
option prices at any time t € [0,T] can be expressed as discounted expected option 
payoffs, where the expectations are computed with respect to a density ft of St- 
The density ft is often called the forward density of St- This paper addresses the 
modeling of the initial density /o and the evolution of the density ft over time, 
{ft}te[o,T]- 

The model is constructed on a filtered probability space {i},J^,{J^t}telo,T],P) 
with expectation operator E. For each t € [0,T], the forward price of a derivative 
pay oS g {St) is the expected payoff computed with respect to the density ft'- 

(1) ngiST) iJ't]^ 1 9{x)ft{x)dx. 



In statements such as the above, to avoid technicalities, the functions mentioned 
are assumed to satisfy measurability and integrability conditions necessary for the 
statements to be meaningful. The market is assumed to consist of n -I- 1 forward 
contracts on European call options with payoffs {St ~ ^"0)+ f^'" *-* ~ 
• • • < Kn- If the original market consists of a mix of European puts and calls, then 
the put-call parity may be used to define an equivalent market consisting entirely 
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of forward contracts on call option payoffs. From ([T]) it follows that the forward 
price processes {Gj}tg[o.T] a-re martingales satisfying the initial condition 

(2) Gi = Jix-Kj) + foix)dx for 0,1,..., n. 

A parametric form for /o will be selected that allows its parameters to be set in a 
straightforward manner from the n + 1 equations in ^ and internally consistent 
forward prices Gg, Gj, . . . , Gg . 

The filtration {J^t}tG[o,T] is assumed to be generated by a standard {n + 1)- 
dimensional Brownian motion {V^, V^, . . . , and we take, for all t, ft to be 

a function with parameters t, V^^ ,Vf^ , . . . , V^^^ that vary over time and other pa- 
rameters that are set in the initial calibration of /g to the current price data. The 
choice of /t allows the R""'"^-valued forward price process (G", G^, . . . , G") to be 
expressed in terms of the Brownian motion (y^, . . , y+i) as 

(G°,Gl,...,G^)-/^t(F,^y,^...,^,"+^) 

for functions ht : R"+^ R"+^, t G [0,r]. It is desirable that the functions ht 
are locally invertible so that the filtration {^/t}tg[o,T] generated by the prices, the 
filtration with an economic interpretation, equals the Brownian filtration. For the 
model to be relevant the functions ht must give rise to price processes with joint 
dynamics that are in line with empirically observed stylized facts for option price 
processes. Moreover, the range of option prices that the model can produce must 
be large enough to capture the fluctuations of observed option prices and avoid the 
need for frequent recalibration. Frequent recalibration of a model's parameters is 
unattractive from a theoretical point of view and limits its practical utility. 

The model for {/t}tG[o,T] set up at time is intended to be relevant also at 
time t > 0. Therefore, it makes sense to require that the realized forward prices at 
time i > should be possible realizations of the model prices Gj , G^ , . . . , G". The 
following example illustrates that a simple model such as Black's model does not 
satisfy this requirement. 

Example 1 (Black's model). Consider the case n = 1 (one forward contract 
and one call option on St) and let {Wt)te[o,T] be standard Brownian motion with 
respect to P. Black's model, see [1], says that 

5T = G°exp{(7oWT- yT} 

which implies that the forward price Gj for the call option payoff {St — K)-^- is 
given by 

Gl = G[]$(di) - i^$(d2), 

The parameter dg solving this equation is the option's implied volatility (implied 
from Gg and Gq). Writing 

2 2 

St = G[! exp {aoWt - ^t} exp {ag(WT -Wt)-^{T- t)} 

and Tt = cr({VFs}s£[o,t]) we notice that the model allows stochastic fluctuations in 
the forward price of St • 

G? = G^xp{aoM^,-^^}. 
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However, the option's implied volatility is required to stay constant over time. In 
particular, the future realized prices are practically guaranteed to violate the model 
which therefore has to be frequently recalibrated to fit the price data. 

One reason for the inability of the dynamic version of Black's model in Example 
[T]to generate future option prices is that the filtration {J-t}t(£[Q,T] is generated by a 
one-dimensional Brownian motion. After the initial calibration the range of possible 
forward prices that the model produces is very limited: it is likely that, after a short 
period of time, the observed option prices lie outside the range of the model. Similar 
problems occur for instance for the local volatility model by Dupire [9] and for many 
stochastic volatility models. In addition, the initial calibration for these models is 
non-trivial. 

In the model we will consider below we want {J't}te[o.T] to be equivalent to the 
filtration generated by the price processes and consider the situation when no price 
process can be determined from the other price processes. 

We do not consider the spot price process for the underlying asset, only its 
value at time T and forward and other derivative contracts written on that value. 
If the asset is a non-dividend paying stock, then the spot price must equal the 
discounted forward price in order to rule out arbitrage opportunities. 

We do not pay attention to the subjective probability views of market partici- 
pants. Therefore it does not make much sense here to discuss equivalent martingale 
measures. However, by requiring that the conditional density process is a martin- 
gale and that it produces realistic dynamics for the price processes we are implicitly 
saying that the model could be a natural candidate for an equivalent martingale 
measure for informed market participants. 

The paper [17j has a similar objective as ours. However, whereas in |17j the 
authors set up a system of stochastic differential equations (diffusion processes) for 
the evolution of the spot price and the implied volatilities and address the diffi- 
cult mathematical problem of determining conditions for the absence of arbitrage 
opportunities, we consider a more explicit but less general class of models for the 
conditional density process {ft}telo,T]- A more general problem is investigated in 
[5] and [13j . where characterizations are provided of arbitrage- free dynamics for 
markets with call options available for all strikes and all maturities. Conditional 
density models, which are studied in this paper, are also studied in |10j . where 
the authors characterize the "volatility processes" {(yl i^)}te[o.T] in the stochastic 
exponential representation 



that generate proper conditional density processes. In contrast, we take a particular 
model for {/t}tG[o,T] as the starting point whereas in [lOj the conditional density 
model is implied from the model for {a{ {x)}t^[o,T]- 

In [4], [8], [Tjj and |14j the authors consider a setting with a finite number of 
traded options for a finite set of maturities on one underlying asset, and characterize 
absence of arbitrage in this setting. Both static arbitrage and arbitrage when 
dynamic trading in the options is allowed are considered. In [4], [7j, and |14j 
explicit Markov martingales are constructed that give perfect initial calibration to 
the observed option prices. 

The outline of this paper is as follows. In Section [5] we consider a rather naive 
model for /o, a distribution of St that reproduces the given option prices, and 
present a straightforward calibration procedure for the model parameters. The 
model is the starting point for the conditional density model for {ft}te[o,T] that is 
presented in Section [3) The theoretical properties of the model and a discussion on 
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how the model can be set up to meet the natural requirements for a good derivative 
pricing model are also included in Section [31 Section S] contains further theoretical 
and numerical investigations of the properties of the conditional density model and 
it is evaluated on S&P 500 index option data and through simulation studies. 

Our contributions can be summarized as follows. We propose a simple model 
for the evolution of the forward density. At each time the forward density is a 
mixture of lognormal distributions which makes it easy to make the initial cali- 
bration of its parameters and price European type derivatives. On a market with 
n liquidly traded call options and a forward contract, the model is driven by an 
{n + l)-dimensional Brownian motion, making it flexible enough to capture realized 
option price fluctuations in a satisfactory way and avoids the need for frequent re- 
calibration. The model is set up so that the filtration generated by the n + 1 price 
processes is essentially, see Section IXTI for details, equal to the (n -I- l)-dimensional 
Brownian filtration. Moreover, the model can easily be set up to capture stylized 
features of option prices, such as a negative correlation between changes in the for- 
ward price and changes in implied volatility. A simulation study and an empirical 
study of call options on the S&P 500 index illustrates that the model provides a 
good fit to option data. 



2. The spot price at maturity 



We start by investigating a very simple model for St, which will be refined 
later, that reproduces the n + 1 observed forward prices. The random variable St 
is assumed to be discrete and takes one of the values < xi < • • • < Xn+2 < oo. 
Let Pq be the forward probability of the event {St = x^}- The initial calibration 
requires solving a linear system of equations of the form Ap ~ b, where p is the 
vector of forward probabilities of the events {St = Xk}'- 

(3) 

/I 1 ... 1 \ 

Xl X2 ... Xn+2 

{Xi - Ki)+ {X2 - Ki)+ . . . {Xn+2 - Ki) + 



( pI \ 

Pa 



/I \ 



Kn)+ J 



V p^' J 



Gg 



V J 



\ {xi - Kn)+ {X2 - Kn)+ ■ ■ ■ (a;,i+2 

If further X2 < Ki, Xn+2 > Kn, and Xk e {Kk~2,Kk^i] for k = 3,...,n-|- 1, 
then the matrix on the left-hand side in ^ is one row operation away from an 
invertible triangular matrix. In particular, the matrix equation Ap — b can be 
solved explicitly for p by backward substitution and then it only remains to verify 
that p is a probability vector. In order to ensure the existence of a probability 
vector solving ([3]) it must be assumed that 

1 _ 

e [0,1] 



(4) 

where we set Kq — 0, and 



- Ctq 



3 



> 1, 



(5) 



K 



Gq 



G: 



>0, j>l. 



The conditions ^ and ([S]) were considered in '6; and ensure that the market of 
linear combinations of forward contracts together with a linear pricing rule is free 
of static arbitrage opportunities. The following result, which is proved at the end 
of the paper, is used as a starting point in the initial calibration of the model for 
the forward price processes presented in Section [31 The result gives (necessary 
and) sufficient conditions for the existence of a discrete distribution of St that is 
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consistent with the forward prices on St- The statement of Proposition [T] below is 
a slight generalization of Proposition 3.1 in 4 . 

Proposition 1 . Suppose that the non-negative forward prices Gq and Gq, . . . , Gg 
on the values St and {St — Kj)^, for j — 1, . . . ,n, at time T > 0, are ordered so 
that Ki < • • • < Kn and satisfy (|4|) and ([5]) . If Xk = Kk-i for fc = 2, . . . , n + 1, 

... . G'',{K2-K,) + GlK,-GlK, 

{K,-K,)-iGl-Gl) ' "'^'^ 

y') ^"+2 ^ t;;^ — — , 

then there exist a unique probability vector (pi, . . . ,Pn+2) such that 

n+2 n+2 

= X! P^^^ '^o = X! Pk{Kk-i - Kj) for j ^ I, . . . , n. 

k=l k=j+2 

The pkS are given by 

Ki+Gl- Gl 



Pi = 

P2 = 

(8) Pk = 



Ki - xi 

xi[Gl -Gl- {K2 -Ki)]+ G"„(K2 -Ki)- GlK^ + Ggi^i 
[Ki - xi){K2 - K^) 

G''^-\Kk-Kk-2) , Gg 



Kk-i~Kk-2 {Kk-i ~ Kk-2){Kk - Kk-i) Kk-Kk-i 
for fc = 3, . . . , n, 

Gg-' Gg(x„+2-i^»-i) 

Kn - Kn-l {Kn ~ Kn-l){Xn+2 - Kn) ' 

Gn 
_ 

Pn+2 — -r^- 

Xn+2 - Kn 

Remark 1. Notice that Q and ^ imply that xi < Ki and that ([7]) implies 
that Xn+2 > Kn ■ Notice also that Proposition [7] says that there exist indicators 
Ik G {0, 1} satisfying Ii + ■ ■ ■ + In+2 = 1 o,nd pk — ^[Ik] such that 



GO = E[^/fea:fc] and G^q = E [( ^ /feXfe - if, 



for j = 1, 



The conditions ([I])-© o-fe sharp: it can be seen from the proof that if any of them 
is violated, then the conclusion of Proposition [I] does not hold. 

Although the model in Proposition [T] for St under the forward probability 
provides explicit expressions for the model parameters in terms of the prices and 
reproduces any set of observed prices satisfying ^ and ([5]) it is not a good model. 
If we want to use the model for pricing new derivative contracts, then wc should 
feel uncomfortable with having a finite grid of points as the only possible values 
for St- For instance, the contract that pays 1 if St takes a value other than one 
of the grid points would be assigned a zero price and this would be viewed as an 
arbitrage opportunity by most (all) market participants. 

A simple extension is to model St as the random variable 

n+2 

(9) 5't = hxkZk 

k=l 

which corresponds to replacing the fixed values Xi, . . . ,Xn+2 by random values 
Xi 'Zi\ , . . . , Xn+2 Zn+2 for some suitably chosen random variables Zi , . . . , Zn+2 that 
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are independent of /i , . . . , In+2 ■ Take 

2 

(10) Zfc = exp{-^T + afeBT}, 

where St is N(0, T)-distributed. ThenE[{xkZk-Kj)+] = G^{xk,(7k,Kj,T), where 

(11) G^(x,CT,ii',T) = a;$(di) - X$(d2), 

rfi — \ — , d2=di-(7VT, 

a\JT ^ 

is Black's formula for the forward price of a European call option maturing at time 
T, where x is the forward price of St, cr is the volatility, and K is the strike price. 

The initial calibration problem for the modified model amounts to finding a 
probability vector p solving the linear equation Ap = b, where A is a square matrix 
with n + 2 rows and columns with 

Ai,k = 1 for all k and Aj_k — G^{xk,<yk, Kj^2,T) for j >2 and all k, 

and where b — (1, Gq, . . . , Gp )"^. The solution p to Ap — b can, as before, be 
expressed as p = A~^b as long as we specify the x^s and akS so that A is invertible. 
In general A will not be close to a diagonal matrix and therefore p = A~^b has to 
be computed numerically. Notice that for a vector b of internally consistent forward 
prices and an invertible matrix A we may find that A~^b has negative components. 
In that case the price vector b is outside the range of price vectors that the model can 
generate. Fortunately it is not hard to determine the range of forward price vectors 
that the model can produce. The simplex S = {p G R""*"^ : P > 0, 1'^p = 1}, where 
l""" = (1, . . . , 1), is a convex set and a linear transformation A of a convex set is a 
convex set. Moreover, the extreme points of S are mapped to the extreme points of 
AS. Therefore it is sufficient to determine the points bk = Ack for k = 1, . . . , n + 2, 
where is the /cth basis vector in the standard basis for R"+^, and investigate the 
convex hull of . . . , 6„+2}. This is the set of price vectors that the model can 
produce. 

3. The forward price processes 

A choice of the initial forward distribution, ^0(2^) = ^{St < x) and forward 
density /o(a;) — Fq{x) has been proposed implicitly from ^ and pH)) . In this 
section, the evolution of the forward distribution and density will be treated as a 
stochastic process {ft}te[o,T], where Ft{x) — P{St < x \ Ft) and ft{x) = Fl{x). 
The filtration {Ft}te[o,T] is taken to be generated by an (n+l)-dimensional standard 
Brownian motion. The n-dimensional Brownian motion corresponding to the first n 
components is denoted by W {W'^ — V'' ior k = 1, . . . ,n) and is used to model the 
indicators /i , . . . , In+2 , whereas the 1-dimensional Brownian motion corresponding 
to the last component is denoted by _B (i? = 1/"+^) and is used to model the 
variables Zi, . . . , Zn+2 as in (fTU)) . 

The forward price at time i of a derivative contract on St with payoff function 
g is given by 

n+2 

H9{St) I ^t] = E [.g(^ /feXfeZfc) I Ft 
fc=i 

n+2 

= ^ P(/fc - 1 I Ft) E[g{xkZk) I Ft]. 
fc=i 
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We consider a partition {Di, . . . , Dn+2} of R" and set Ik = I{Wt G Dk}- The 
factors P(/fe = 1 I J^t) and Fj[g{xkZk) \ J^t] can be computed as follows: 

Pih = 1 I ^0 = PiWr e Dk I m), 

2 

E[g(a;feZfc) I = E [3(2:^ exp { " + ^feST}) | 

We write 

= P(4 = 1 I ^t) = PiWt + WT-WteDk\ Wt) = , 



where $„ is the standard Gaussian distribution in R". Note that the stochastic 
process {pt}te[Q,T]j where pt — {pj, . . . is a martingale on the simplex S = 

{p G R"+^ ■ P > 0, i^P = 1} with the property that pr G {ei, . . . , e„+2}, where 
the efeS are the basis vectors of the standard Euclidean basis in R"+^ . The forward 
prices at time t are given by 



G'i=J2ptx'^ and =^p^GB(x^a,,if„r-t) for j = 1, 



.,n, 



where denotes Black's formula ([TT|) for the forward price of a European call 
option and 



Xt =XkE 



2 2 

exp { - + akBr} \ St] = Xk exp { " + f^fc^t}- 



3.1. Tracking the Brownian particle in continuous time. In order to use 
the model at time t G (0, T) for pricing a European derivative with payoff function 
g, it is necessary to know the location of the Brownian particle {W^, . . . , W",Bt). 
That is, given the observed forward prices (Gj , . . . , G") we need to infer the lo- 
cation of (Wf, . . . , Wp, Bt). We may express (G^ , Gj , . . . , G") as the value of a 
function ht evaluated at [W^ , . . . , Wp,Bt). The filtration {5t}tg[o,T] generated by 
the vector {G'^,G^, . . . , G") of price processes is therefore smaller than or equal to 
the Brownian filtration {J^t}te[o,T] generated by {W^, . . . , VF", B). We now inves- 
tigate the functions ht in order to compare the two filtrations and to determine the 
dynamics of the price processes. 

The mixture probabilities can be written as p^ — p^{Wt), where 



(12) p^{w)= I (27r(r-t))-"/2exp{ 



[x — w)^ {x — w) ' 
2{T-t) . 



■dx. 



Write ht = (/i°, . . . , /i"). Then the forward price Gj can be expressed as Gj 
h°{Wt,Bt), where 



n+2 



h°{w,b) = ^p^{w)xkexp\^- ^t + akb}. 



k=l 



Similarly, Gl = hi {Wt,Bt), j — 1, . . . ,n, where 



n+2 



hi{w,b) = Y,PtMG^[xkexp{- ^t + akbY<Jk,K,,T~t 

k=l 

n+2 2 

^pj^(w)(xfcexp{ - + akb^^idi) - K^<S?{d2) 



k=l 
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with di = di (j, k,b,T — t) and d2 ~ ^2 (j, k,b,T ~ t) given by 
di = 7=^( - + <^kh + \og{xk/K,)) + 



d2 = di - ak^/T - t. 
In particular, 

(13) {GlGl...,G^) = ht{Wt,Bt) 

= {h'l{Wt,Bt),h\(Wt,Bt), hnWt,Bt)). 

If, for every t G [0,T), ht : R""*"^ — R"+^ is locally one-to-one everywhere, then 
an (n -I- l)-diniensional trajectory for the forward prices can be transformed into a 
unique (n-l-l)-dimensional trajectory for the (n+l)-diniensional standard Brownian 
motion (W, B). From the inverse function theorem (Theorem 9.24 in fl6|) we know 
that if the Jacobian matrix 

d^i^^b) ... g^{w,b) ^{w,b) > 

(14) K{w,b)^ : : 

of the continuously differentiable function ht is invertible at the point {w,b), then 
ht is one-to-one in a neighborhood of (w, b) and has a continuously differentiable 
inverse in a neighborhood of ht{w, b). The set 

Tt = {{w,b) £ R"+i : det h[{w,b) =0} 

is the subset of R"+^ where ht is not locally one-to-one. 

In order to investigate the sets Tt and in order to express the dynamics of the 
price processes using Ito's formula the partial derivatives of the functions ht must 
be computed. We find that 

-{w,b) = ^Pt{w)akXkexpi^ - -^t + '^kbj 



814. 
db 

k=l 



and 



dh' "+2 , ^2 . 

-griw,b) = ^Pt^(u;)$(di)crfcXfcexp| - -^t + akbj, 

k=l 

where di — di(j, k) depends on j and k through Kj and crfc. Similarly, 

and 

Finally, 



Numerical investigations, illustrated in Figure[TJ indicate that Ft is a smooth surface 
of dimension n that varies continuously with t. If the function det h't : R"+^ — )■ R 
has a nonzero gradient almost everywhere in Fj = {(w', ^) G R""''"'^ : det/ij = 0}, 
then the implicit function theorem (Theorem 9.28 in |16j ) implies that Ft is a con- 
tinuously differentiable hypersurface in R"+^. Similarly, if the function {t,w,b) t-^ 
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det h^{w,b) has a nonzero gradient almost everywhere in F = {{t,w,b) E R"+^ : 
det h'^ = 0}, then F is a continuously differentiable hypersurface in R"+^ from which 
we conclude that the Fjs vary continuously with t. If the gradients are nonzero al- 
most everywhere, then we conclude that P{{Wt, Bt) S Fj) = for all t but that 
P{{Wt,Bt) e Fi for some t) > 0. In particular, if 

T = M{t>0:{Wt,Bt)ert}, 

the first time that the (ri+l)-dimensional Brownian motion (W, B) arrives at a point 
where h^- is not locally invertible, then the trajectory of {(VFt, S()}ig[Q ^.j is uniquely 
determined by the trajectory of {hf {Wt, Bt)}t£io,T]- Therefore, r is a stopping time 
with respect to {Gf}t£io,T] and GtAr — J^tAr- However, whether the trajectory of 
{{Wt,Bt)}t(z[o^T] is uniquely determined by the trajectory of {ht{Wt, Bt)}te[o,T] or 
not depends on the function Ht in a neighborhood of {Wr,Br) G Tt- 

In practice, only discrete observations of the forward prices are available, so it 
will be impossible to track the Brownian motion exactly based on the discretely 
observed forward prices. This issue is treated in some detail in Section S] where 
both a local linear approximation and a particle filtering method is applied to track 
the location of the Brownian particle. 




Figure 1. Contour plots of the Jacobian determinant det/ij, for 
n = 2 and 6 = 0, as a function of (wi, W2) at times t = 0, t = 0.5, 
and t — 0.9. The zeros of the determinant are displayed along the 
dotted curves. The functions ht correspond to a forward density 
process calibrated to S&P 500 option data presented in Section S] 
and parameterized as in (fT5)) . 

3.2. The forvirard price dynamics. Many popular models for derivative 
pricing are based on modeling the dynamics of the underlying spot price or forward 
price directly. Examples are Black's model, Dupire's model, and stochastic volatil- 
ity models. Our starting point is a model for the dynamics of the forward density. 
From the model for the forward density process, the dynamics of the forward price 
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process {Gj}tg[o.T] a-re derived from the expressions for the partial derivatives of 
ht and Ito's formula (Theorem 33, p. 81, in 15 ): 



n n+2 pt ft k 2 

= Gg + EE/ /^(W^.Kexp{-^. + afeS,}dW^: 

1=1 fc=l "'0 
n+2 2 

+ H / Ps(W^s)f^fea;feexp| - ^s + crfcBs}dS^ 

A;=l 2=1 * 



From e.g. the martingale representation theorem (Theorem 43, p. 186, in 115' ) it 
follows that the last sum above vanishes so that 

n n+2 nt n k 2 

= + ^ ^ / f^iWs)xk exp {~^s + akB,}dW: 

1 = 1 fe=l "'0 
n+2 „t 2 

+ H / Ps(W^^*)o-fca;fcexp| - ^s + crfcS^jdS^. 
fc=l "^0 



The derivatives computed so far can also be used to study the conditional 
density process {ftix)}te[o,T]- The conditional density 



n+2 



fc=l 



is a convex combination, with random probability weights as above, of lognormal 
densities ft{x), where 

Ito's formula and the martingale representation theorem yield, where the depen- 
dence of ft{x) on Wt through the pj's and on Bt through the f^s has been sup- 
pressed, 



_^ /o 



dwf ' ds 
" "+2 /•* ^i^fc 



i=l fc=l "'O ' 

/^;/.n.)(^2SWfi)±££Zi^),B. 

^ c7k(T-s) ) 
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Since ft {x) > everywhere we may write ft (x) as a stochastic exponential 



(15) 

where (Vi, 



ft{x) = fo{x) + / {<y!nx)f.{x)dv:, 

{Wi,...,Wn,B) and 
{af)\x)^fs{xy 



,Vr 



n+1) 



{air+\x) ^ f^ix) 



k=l 



n+2 
fc=l 



\og{x/xk)+alT/2-akB, 
<Jk{T-s) 



Conditional density models of the form (ITSI) are studied in |10j . In |10l , the authors 
characterize the processes a{ {x) that imply that {ft{x)} in (fT5)) is a conditional 
density process (a forward density process with the choice of numeraire considered 
here), and provide several explicit examples. Here, we consider a particular forward 
density process and determine the corresponding volatility process {a{ (x)}. 

3.3. Explicit computations in the case n = 2. The expression for p^, 
k = 1, . . . ,n + 2 in ((T^ is an integral of a Gaussian density over a set Dk- In this 
section explicit evaluation of the partial derivatives of will be performed in the 
case where n = 2 and Dk is a cone. 

Recall that the case n = 2 corresponds to one forward contract on St and 
two call options on St- In this case W = {Wx,Wy) is a Brownian motion in 
R^. We choose the sets Di, D2, D^, D4 as cones or unions of cones because that 
gives a convenient parameterization for numerical computations and because the 
configuration of the number of cones and their placement can be rather easily 
modified to produce dynamics for the price processes that we find reasonable. 

Let D be the cone in the first quadrant expressed in polar coordinates as {(r, v) : 
r > 0,v G [0,6']} for 6 S [0,7r/2]. With Z denoting a random vector with the 
standard two-dimensional Normal distribution we write 



(16) 



p{w + Vt^z e D) 
1 



D MT ~ t) 



exp 



1{X- + {y- Wyf 



T-t 



^dxdy 



1 



■ exp 



The derivative of (H 



y/tane 2'k{T - t) dw^ 

1 X 



exp 



with respect to Wx is 

1 d r 1 (x - 

I 2" 



\{x- + (y - WyY 



T-t 



'Jx? + {y-wyf 



^dxdy. 



T-t 



^dxdy 



1 



27r(T - t) 

^''^ = 2^) 
The identity 



y/,,^eMT-t) T-t 

1 {x 
2~ 



exp 



{ 



l{x- Wxf + {y- WyY 



T-t 



^dxdy 



exp 



exp 



{ 



T-t 

1 {y/i&nO-Wxf + 



y / tan f 



dy 



(y- 



T-t 



"-^}dy. 



{ay-b)' + {y-cf 



(y _ O-b+A 



{ac-bf 
l + a2 
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with a = 1/ tan^, h — Wx and c — Wy can be used to write the integral in (jl7p as 

exP| - 2(T-t)(l + l/tan-e)/ \ 2 ^^^^^ / 



V2.(r-t)(l + l/tan^^?) 7o V'2vr,^^ 



The integral expression may not look pretty but can be written explicitly as 

Q , 

P{w + VT-tZ eD) 



dwr, 



( _ 1 (Miy/ tan \ 

'\ 2(T-t)(i+i/taii2 0)j / Wx/ tan9 + Wy 



exp 



v/27r(T-t)(l + l/tan2 6*) ^ 0(1 + 1/tan 

in terms of the univariate standard Normal distribution function $. Similar com- 
putations for the derivative of (jl6p with respect to Wy give 

^ p(.« + VT^z e D) = !!4=sici>(^) 

^ ^2Tr{T - 1) \VT^J 

~ 2 (T-t)(i+tan2 9) j ^ / + ti;„ tan ( 



/ Wx + Wy tan t/ \ 



V27r(T-<)(l +tan2 0) ^ ^(y _ + tan^ 6*) 

Let Z? be a cone that can be expressed, in polar coordinates, as {(r, w) : r > 0, u € 
[0,0 + 9]}, where 6 € (0, 7r/2] and (/> S [0,27r — 9]. Let further be the matrix 
corresponding to a clock-wise rotation of angle </> around the origin so that O^D is 
of the form, in polar coordinates, {(r, w) : r > 0, u G [0, 9]}. Then 



P{w + VT-tZ eD)= P{0,pw + VT-tZ e O^D) 

and the above computation, with w — O^w instead of w, can be used to compute 
the partial derivatives of P{w + \JT — tZ e D) with respect to Wx and Wy. With 
u = 0(f,ei and v = 0^62, where ei and 62 are the standard basis vectors in R^, we 
get 



P(u; + VT^Z e D) = cos (/)—- P{w + y/T-tZ G O^D) 



dwx dw 



d 



sin 0- — P{w + y/T -tZ e Oa,D) 

OWy 



and similarly 



P{w + VT-tZ eD)= sin 0--^ P{w + y/T-tZ e O^D) 



dwy dwx 



d 



+ cos (f)—^ P{w + y/T-tZ e OaD). 

OWy 

4. Calibration and evaluation of the model 

To calibrate and evaluate the model we use 41 daily adjusted closing prices over 
a 59 day period from September 22nd 2011 to November 19th 2011 of European 
put and call options with strike prices 1150, 1175, 1200, 1225, and 1250 on the S&P 
500 index value. The options mature on November 19th 2011. A rather short 
time series of option prices is selected in order to have price data corresponding to 
sufficiently large traded volumes so that the option prices can be considered to be 
relevant market prices at the end of each trading day. The risk-free interest rate 
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is set to 0.5% (corresponding approximately to the three-month LIBOR rate) and 
the put-call parity 

Ct(i^)-Pt(AO = e-"°"5'^-*'(GS'-if), A' e {1150,..., 1250}, 

is used, for the most traded pair of put and call options on each trading day, to 
calculate the forward prices Gj of the underlying asset. For example, on September 
22nd 2011, the put and call options with strike 1150 were the most traded options 
and their option prices were used to calculate the initial forward prices for delivery 
of the value of the S&P 500 index on November 19th. The initial forward price 
was calculated to Gq = 1128.12. During the analyzed time period, September 22 
- November 19, the forward price increased. Simultaneously the largest trading 
volumes shifted from the options with strike 1150 to the options with strike 1200. 

4.1. The initial calibration and model specification. In this section the 
initial calibration and model specification for the S&P 500 options will be explained 
in some detail. The first step is to select the grid parameters Xk, the volatilities Ufe, 
and the partitions Dk, k = 1, . . . , n -|- 2. The parameters will be selected to get a 
reasonable shape of the initial density /o and such that the evolution of the prices 
have features that are present in real data. 

Let us start by considering the initial density /q. As a reference density we 
will consider the density, , resulting from Black's formula with a fitted volatility 
smile. At time (September 22nd 2011) the implied volatilities, af , i = 1, . . . ,n, 
corresponding to the strikes Ki, . . . , Ar„, are computed using Black's formula for 
European call option prices, 

C^{K) = e-o ™5?'(G°<i>(di) - K^d2)), 

log(G°/X) a{K)VT /- 

A second-degree polynomial (volatility smile or volatility skew) is fitted to the 
implied volatilities and the formula 

d'^C^{K) _ -0.005-59/365 

dK^ 90 (^), 

see e.g. [2] or gives the probability density for St implied by the volatility 
smile and Black's call option price formula (see e.g. |18| for details). The prices of 
the call options with the strike prices 1150, 1175, 1200, 1225, and 1250 produce the 
volatility smile (second degree polynomial) shown in the upper left plot in Figure 
[5] Notice that the fit is rather poor. The corresponding implied density is 
shown in the upper right plot in Figure [2] A closer look at the data reveals that 
the call option with strike price 1175 only has 82 registered trades, so the price 
of that contract may be unreliable. If that implied volatility is omitted, then the 
volatility smile in the lower left plot in Figure [5] and the implied probability density 
in the lower right plot in Figure [5] are obtained. Note that the resulting probability 
density is smooth, unimodal and left-skewed. We will assume that the volatility 
smile in the lower left plot in Figure the graph of the second-degree polynomial 
K I— > a^{K) fitted to the implied volatilities, corresponds to correct market prices 
which will be used in the calibration of the model in Sections \4.1.1\ and \4.1.2\ 

As can be seen in the lower right plot in Figure [U implied volatilities that are 
decreasing with the strike price correspond to a left-skewed implied density function 
for St- This observation is in line with much of the empirical analysis on option 
price data, see e.g. 

To select the parameters in our model for /o, the parameters ak will be chosen 
to produce a left-skewed implied density. A commonly held view is that the changes 



14 



H. HULT, F. LINDSKOG, AND J. NYKVIST 



in the implied Black's model volatility and log-returns of the forward (or spot) price 
are negatively correlated, corresponding to different market responses to good and 
bad stock market information, see e.g. |12) . This behavior is in line with our 
findings based on the small option price sample used here: for the options with 
strike prices Ki = 1150 and K2 = 1200 the sample correlations between the daily 
log-returns \og{G^j^i / G^) and the implied volatility changes (tP(_|_;^ — a^^, i = 1,2, 
are both —0.35. For this reason it makes sense to choose the parameters at to be 
decreasing in k. If the probability mass of the p^s is shifted towards lower indices 
k, then the forward price decreases and the implied volatility increases. 

In Sections 14. 1 . II and 14. 1 .21 below we present the initial calibration of the model 
to n = 2 and n = 5 option prices, respectively. In both cases it is assumed that 
the volatility smile in the lower left plot in Figure [5] corresponds to correct market 
prices. 

4.1.1. Initial calibration with two options, n — 2. Let us first consider two 
options, n = 2, with strikes Ki = 1150 and K2 — 1200. The corresponding 
implied volatilities are cr^(ifi) = 0.33082 and a^{K2) = 0.29777. Since the current 
forward price of the underlying asset and the risk-free interest rate are known 
we can use Black's formula to calculate the option prices. We get Cq = 49.575 
and Cq = 26.434, and the corresponding forward prices Gq = 49.615 and Gg = 
26.455, respectively, of the options. To calibrate the model to these prices, the 
parameters xi, X4 and ak, k = 1,...,4 need to be specified so that the vector 
Po = (Poj • • • iPo)'^ is ^ probability vector, i.e. has non- negative components that 
sum up to one. Proposition[l]is used to obtain x™^^ — 1016.81 and a;™'" — 1257.11. 
We choose xi = 950 and X4 = 1300 and initially choose at = 0.01 for all k. The 
resulting density function fa{x) — /o(a;; cti, . . . , CT4) is displayed in the upper left 
plot in Figure [3] Even though these parameters are consistent with the observed 
prices, we are not comfortable with the appearance of the resulting probability 
density for St- To get a smoother density we need to increase the CfeS. First we 
increase the o-feS simultaneously as long as po stays a probability vector. It turns 
out that fj/c = 0.0542 is the largest possible value, but the corresponding density 
fo{x) = fo{x; tJi, . . . , (74) is not left-skewed. Since the lognormal density is right- 
skewed, the natural approach is to increase the cTfeS for small fcs and decrease the 
CTfcS for large ks. The parameters ci, 0-4 = 0.18, 0.08, 0.06, 0.03 give the density 
fo{x) = /o(x; CTi , . . . , CT4) in the upper right plot in Figure[31 which is rather similar 
to the implied density . We summaries the chosen model parameters: 

( xi = 950, X2 = Ki = 1150, X3 = K2 = 1200, X4 = 1300, 
(18) 71 = 2: <^ cri,...,cr4 = 0.18,0.08,0.06,0.03, 
[ pI,...,P^ w 0.29,0.14,0.51,0.07. 

4.1.2. Initial calibration with five options, n — 5. Here n = 5 options, with 
strikes 1100, 1150, 1200, 1250, 1300, are considered to illustrate that the calibration 
procedure easily handles more than two option contracts. Similar to the setting in 
Section 14.1.2 1 the parameters xi, xr and ak, k — 1, ... ,7 are specified so that po is a 
probability vector. Using Proposition[T]to obtain x™^^ — 968.86 and 2;™'" — 1321.8. 
We choose xi = 950 and xi — 1400 and begin by choosing ak = 0.01 for all 
k. The resulting probability density is displayed in the lower left plot in Figure 
[31 We would like the density to spread out the probability mass more evenly 
and therefore we increase the akS simultaneously as long as po stays a probability 
vector. It turns out that ak — 0.027685 is the maximum possible value. Increasing 
the CTfcS for small A:s and decreasing the akS for large fcs will produce a left-skewed 
density. The density shown in the lower right plot in Figure [3] corresponds to 
ai,...,ar = 0.21,0.045,0.028,0.025,0.025,0.02,0.01. We summaries the chosen 
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Figure 2. The upper left plot shows the miplied volatilities for 
the five call options, o, and the fitted second-degree polynomial 
(solid curve). The upper right plot shows the probability density 
for St derived from the volatility smile in the upper left plot. 
The lower left plot shows the implied volatilities for the four call 
options, o, and the fitted second-degree polynomial (solid curve). 
The lower right plot shows the probability density for St derived 
from the volatility smile in the lower left plot. 

model parameters: 
(19) 

{xi = 950, X2 =Ki = 1100, X3^K2 = 1150, X4 = K3 = 1200, 
X5 = = 1250, X6 = K5 = 1300, X7 = 1400, 
CTi, . . . , 0-7 = 0.21, 0.045, 0.028, 0.025, 0.025, 0.02, 0.01, 
pl,...,plRi 0.26, 0.23, 0.08, 0.15, 0.15, 0.13, 0.002. 

In principle the model can be set up and calibrated to an arbitrarily large number 
of option contracts. In practice, however, it is difficult to find a large number of 
reliable option prices for a wide range of strikes. For example, we notice that in 
our data the options that are actively traded all have strikes close to the current 
spot price of the underlying asset. Options with strike prices that are far from the 
current spot price have none or very few trades, so their daily closing prices are 
unreliable. 

4.1.3. Model dynamics and selection of sets Dk- To examine the dynamics of 
the model we consider the model parameterized as in ^TE\\ . The sets Di, . . . , D4 are 
chosen as cones and placed in increasing order starting at the x-axis. The cones 
are illustrated in Figure [TJ Then, N = 5000 trajectories are simulated of the 3- 
dimensional Brownian motion (Gaussian random walk) {W, B) with 50 time steps 
for each trajectory, corresponding roughly to the number of days of the sample of 
option prices. For each trajectory and each time step the corresponding forward 
price of the underlying asset and the forward prices of the two call options with 
strike prices Ki — 1150 and K2 = 1200, respectively, are calculated. Next, for 
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Figure 3. The upper plots show the model density fo{x) = 
foix; CTi, . . . , (74) (solid curves) and the implied Black's model den- 
sity Qoix) (dashed curves). The upper left plot corresponds to 
CTfc = 0.01 for all k and the upper right plot corresponds to 
ai,...,ai = 0.18, 0.08, 0.06, 0.03. The lower plots show the model 
density /o(a;) = /o(a^; fi, . . . , 0-7) (solid curves) and the implied 
Black's model density q^{x) (dashed curves). The lower left plot 
corresponds to ak — 0.01 for all k and the lower right plot corre- 
sponds to CTi, . . . , (77 = 0.21, 0.045, 0.028, 0.025, 0.025, 0.02, 0.01. 

each trajectory and each time step Black's formula is used to calculate the implied 
volatilities for the two options. Finally, the correlation between price changes of 
the underlying asset and changes in the implied volatilities is calculated. The 
histograms in Figure S] reveal that the correlation is negative in most simulations 
with mean values —0.51 and —0.56. These values appear to be in line with empirical 
studies, e.g. |12] . 

4.2. Tracking the Brownian particle in discrete time. As observed in 
Section [XT] the vector of forward prices at every time t G [0, T] can be expressed as 
(G?, Gi, . . . , GJ')^ = htiWt^,..., Wp, Bt), where W^,...,W'',B are independent 
one-dimensional Brownian motions. If the vectors yt of observed prices are within 
the range of the model, then yt — ht{xt) for all t E [0,T], where the function ht 
is given by and {xt}ti£[o,T] is an observation of the trajectory of the (n -|- 1)- 
dimensional Brownian motion. Recall that ht is locally invertible for t smaller than 
a stopping time r measurable with respect to the filtration generated by the price 
processes. In principle it is possible to uniquely determine the trajectory {xt}t£[o,T] 
from that of the price process {yt}te[o,T]- 

In practice, the situation is more complicated because the price data consist of 
daily closing prices {ykAt : k — 1,. . . ,T/At}. In particular, the local one-to-one 
property of the functions ht do not guarantee that the trajectories of the driving 
Brownian motions can be well estimated. In this section the aim is to estimate 
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Figure 4. The plots show histograms for the sample correlations 
between log-returns of the forward price of the underlying asset 
and the implied volatility changes, based on 5000 simulated price 
trajectories. The left plot corresponds to implied volatilities for 
the call option with strike price 1150 and the right plot for the 
strike price 1200. 

the Brownian motion at the observation times, {xkAt '■ k — 1, . . . , T/At}, from the 
observed option prices. 

4.2.1. Local linear approximations. Since the function ht in (fT^ is continuously 
differentiable we may approximate ht{x) in a neighborhood of a point Xq by the 
best linear approximation ht{xo) + h'^{xo){x — xq)- Since the htS are, up to time 
T, locally invertible we may use the linear approximations of the hfS together with 
observations ykAt of the prices to obtain, iteratively, estimates XkAt of the Gaussian 
random walk XkAt'- 

(20) J(fe+l)At = XkAt + [h[k+l)AtixkAt)]^^ {y(k+l)At - VkAt), Xq = 0. 

However, the time step Ai corresponding to daily prices is rather large which implies 
that the linear approximation may be inaccurate. Moreover, the Jacobian matrices 
^(fc+i)At(^fe'^t) be too close to singular leading to poor estimates of the x^At^- 

4.2.2. Particle filtering. An alternative approach to the local linear approxi- 
mation is to use an auxiliary particle filter to estimate the XfeAtS, or rather the 
posterior distribution of the XkAt&- The particle filtering approach considered here 
works as follows. 

(1) At time kAt we have R particles at locations a\^^, ct^At^ where = 
for all j. 

(2) To each particle a^^j a first-stage weight \j is assigned, given by 



Aj = , where 9^ = (t>{y(k+i)Au ^(fc+i)At(aiAt)' ^i) 

2^,:=i ^» 

and (/)(y; /i, E) denotes the density at y of the Normal distribution with 
mean /i and covariance matrix E. 

(3) Draw with replacement from the index set {1, . . . according to the 
weights Aj to produce R indices ni, ...,n/f. 

(4) For each j set 5^j,^x)At ~ '-'^fcAt ^" V AtZj , where the ZjS are independent 
and standard Normally distributed vectors. 

(5) To each particle oP{^k+i)At ^ second-stage weight tt^ is assigned, given by 

, 0(y(fe-M)At;^(fe+i)Ai(a(fc+i)At)>^2) 

TTj — — -jT^ , where Wj = „ . — . 

E^li Wi 0(y(fe-n)At; /i(fc+i)Ai(afcAt)> ^i) 
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(6) Draw with replacement from the set {c(lk+i)At' ■ ■ ■ ^'^fk+i)At} according 
to the weights iTj to produce the set of particles {a(i.+i)At' ■ • • > '^(fc+i)At}- 
In order to use the particle filter the filter parameters R, Si, and S2 must be spec- 
ified. Notice that for each k the sample . . . , a^^d forms an empirical distri- 
bution that approximates the conditional distribution ofxkAt given yo, yAt, • • • , UkAt- 




so 100 150 200 2B0 300 3S0 400 450 SOO ' 50 100 150 200 250 300 350 400 450 500 



Figure 5. The left column shows the resulting filter distributions 
for the X-, y-, and 6-coordinate respectively. The right column 
shows the resulting estimates from the linearization. 

4.2.3. A simulation study. In this section the performance of the particle filter 
and the local linear approximation will be illustrated in a small simulation study. 

Consider the model for n = 2 with parameters given by ([T5|. Take T = 1 
and At = 1/500 and simulate {(GfeAt, G^ac G'^At) ■ ^ = 1,...,500} by feeding 
the model with a Gaussian random walk whose increment distribution is the 3- 
dimensional Normal distribution N3(0, Ai/), where I denotes the identity matrix. 
The problem we consider here is to estimate the location of the Gaussian random 
walk from the simulated price data {{GkAt, G\^^,G\^^) : k — 1, . . . , 500}. 

The particle filter parameters are selected as follows. The number of particles is 
R — 250. The matrix Ei is chosen as the sample covariance matrix of the simulated 
increments for the 3-dimensional forward price process. The matrix E2 needed to 
assign the particles' second-stage weights is chosen as E2 = Si. 
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Figure 6. The left column shows the true price trajectories 
(black) and the resulting filter prices (gray) for G, Gi and G2 
respectively. The right column shows the corresponding price tra- 
jectories from the linearization. 



The output of the particle filter is a distribution of the location of the three- 
dimensional Gaussian random walk. The three components of the true simulated 
random walk and the corresponding estimates from the particle filter are displayed 
in the left plots in Figure [5] The empirical distributions of the particles estimating 
the location of the random walk are displayed in grey on top of the true simulated 
trajectories. The estimates from the local linear approximation of the random 
walk trajectories are displayed in the right plots in Figure [5l The particle filter 
approach is reasonably good at tracking the underlying Gaussian walk, whereas 
the performance of the local linear approximation is clearly worse. 

In addition, the forward prices of the index and the two options are recalcu- 
lated using the corresponding particle filter estimates and linearization estimates, 
respectively. For the particle filter, at any given time each particle (an estimate 
of the location of the Gaussian random walk) gives rise to a forward price and 
the weighted sum of the prices corresponding to different particles is the value of 
the gray price trajectory in the left plots in Figure [51 The true price trajectory is 
plotted in black (the one corresponding to the simulated Gaussian random walk). 
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The two price trajectories are essentially indistinguishable. The plots to the right 
in Figure El show the price trajectories computed from the linearization estimates 
of the Gaussian random walk (in gray) and the true price trajectories (in black). 
The linearization estimates of the Gaussian random walk do not reproduce the 
simulated price trajectories as accurately as the particle filter estimates. 

4.2.4. Tracking the Brownian particle for S&P 500 option data. Now that the 
initial calibration of the model and the particle filter is well understood, the particle 
filter is applied to the S&P 500 option price data; 4f vectors of daily closing prices 
for the index forward and two call options. The results with n = 2 are shown in 
Figure [71 The particle filter distribution of the underlying Gaussian random walk is 
rather wide but nevertheless are the computed prices based on the filter estimates 
very close to the real prices. The plots in Figure [7] demonstrate that the model is 
very good at reproducing the true price trajectories. 

For each of the times t E {0, 10,20,30} days from today, the set of particles 
from the particle filter is used to compute model prices for a fine grid of strikes. For 
each of these strikes a call option price is computed as a weighted average (second 
stage weights) of the model prices corresponding to different particles. Then, the 
produced prices are transformed into implied volatilities using Black's formula. 
For each of the four times, the procedure thus produces a volatility smile (a set 
of implied volatilities), and we observe how the volatility smile varies over time. 
Figure [5] shows that the model and particle filter produce volatility smiles at all 
times that appear to be reasonable. 



Appendix A. Proof of Proposition [T] 

Consider the equation Ap = & in ([3]) with the choice of XkS according to the 
statement of the proposition. Clearly, it has a unique solution. We need to deter- 
mine this solution p and verify that p G [0, 1]"+^. Using backward substitution, we 
solve for the three last probabilities Pn+2,Pn+i,Pn to obtain 



Xn+2 ~ Kn ' 

Gg-' Gg(.T„+2-i^n-l) 

Kn - Kn-1 {Kn ^ Kn-l)(Xn+2 - Kn) ' 

Gg-^ G'^^-\Kn~Kn-2) ^ Gg 

Kn-l — Kn-2 {Kn-l ~ Kn-2){Kn — Kn~l) Kn — Kn-1 

We begin by showing, by a standard induction argument, that pk can be written 
as in for k — 3, . . . ,n. We know that this holds for fc = n. We now assume that 
it holds for k = n — j,n — j + 1, n and show that it holds for k = n — j — 1. We 
know that 

n+l 

Gq~^~^ = ^ Pk{Kk-l - Kn-j-l) +Pn+2{Xn+2 ~ Kn-j-l) 
k—n—j+1 

which means that 

_ Gq^^ - (Kn-j-l - Kn-j-z)Pn-j Pn+2{Xn+2 - Kn-j-z) 

Pn-j-1 — V 



Pn+2 = 
Pn+1 = 
Pn = 
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Figure 7. The upper plots and the left middle plot show the em- 
pirical distributions (in gray) of the particles of the particle filter 
applied to real price data, and the mean of the empirical distribu- 
tions (in black) for the x-, y- and 5-coordinate, respectively. The 
right middle plot and the lower plots show the real price trajectories 
(in black) for G°, and and the price trajectories computed 
from the filter estimates (in gray). 



Inserting the expressions for Pn+iiPn+2 and the expression for pk for k = n — j, 
j + I, ...,n, and collecting the terms we obtain 



Kn-j-2 — Kn-j-'i {Kn-j-2 — Kn-j-3)iKn-j-l — Kn-j-2) 

_! ^fo 

Kn-j-l — Kn-j-2 
Kn-j-2 — Kn-j-3 Kn-j-2 — Kn-j-3 
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Figure 8. The plots display volatility smiles produced by the par- 
ticle filter and the model fitted to the S&P 500 option data, at 
times t = days (upper left), t = 10 days (upper right), t = 20 
days (lower left), and t = 30 days (lower right). 



where 



7n-fc = 



Kn-k~l — Kn~k~3 _ [Kn-k — J^n-fc-s) (-^n-fc+l — -^K-fc-l) 
Kn-k ~ Kn-k-1 {Kn-k — Kn,-k-l){Kn-k+l — Kn-k) 

Kn-k+1 — Kn-k-3 



Kn-k+1 — Kn-k 

for fc = 1 , . . . , J , and 

_ Kn-l — Kn-j-3 _ {Kn - -^n-j-s) (a^ri+2 — J^n-l) 

^ Xn+2 — Kn-j-3 
Xn+2 - Kn 

Straightforward calculations show that jn-k = and 7„ = 0. By induction we 
have therefore shown that ([8]) holds for fc = 3, . . . , n and it remains to solve for pi 
and p2- We have 

n+2 

Pi + P2 = 1 - ^ Pfc 

k=3 
n+1 

PlXl +P2K1 = Go - ^PkKk-l -pn+2Xn+2- 
fe=3 

Using dS]) we obtain 



Pi + P2 = 1 - 



rii /^2 

K2-K1 



PlXl + P2K1 = Gg 



GIK2 ~ GlKi 



K2~Ki 
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Solving for pi and p2 gives 



Pi = 

P2 



Ki~ xi 

xi[{Gl - Gl) - {K2 - ifi)] + Gl{K2 - Ki) - GIK2 + GlKi 



We must show that (pi, . . . ,Pn+2) corresponds to a probabihty distribution, i.e. that 
Pk >0 for aU k (and pi + ■ ■ ■ + Pn+2 = !)• 

Pi > : Notice that pi > is equivalent to if 1 + Gq — Gq > which follows from 
©■ 

Pk > for k = 3, n : We have shown that pk is given by ([5]) ioi k — 3, . . . ,n, 
i.e. that 

1 f^<k-2 Kk - Kk^2 ^k-i , Kk-i-Kk- 



Pk 



A^o - Kk-Kk-^' + Kk-Kk-^ 



Kk-i - Kk-2 

The non- negativity of pk is therefore an immediate consequence of ([S]) . 
Pn+i,Pn+2 > : Notice that 

^ Gg-^(x„+2 - Kg) ~ Gg(x„+2 - Kn-l) ^ Q 
{Kn - if„-i)(a;„+2 - Kn) ~ 

is equivalent to GQ~^(a;„+2 — Kn) — Gq (a;„+2 — Kn-i) > 0. Solving for Xn+2 shows 
that the latter is equivalent to 



(21) Xn+2 > 



Gq — GqK„ 



Moreover, since Xn+2 > Kn it obviously holds that Pn+2 — Gq / {xn+2 — Kn) > 0. 
P2 > : First note that P2 > is equivalent to 

xiiK2 ^Ki + Gl- Gl) < G°oiK2 - Ki) + GIK^ - GJX2. 

Moreover, from @ we know that K2 — if 1 + Gq — Gj > and therefore p2 > is 
equivalent to 

. Gl{K2-K,) + GlK,-GlK2 
^ ^ K2-K, + Gl-Gl 

□ 
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